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Abstract 

We find the light-cone wavefunction representations of the Sivers and Boer-Mulders distribution 
functions. A necessary condition for the existence of these representations is that the light-cone 
wavefunctions have complex phases. We induce the complex phases by incorporating the final- 
state interactions into the light-cone wavefunctions. For the scalar and axial- vector diquark models 
for nucleon, we calculate explicitly the Sivers and Boer-Mulders distribution functions from the 
light-cone wavefunction representations. We obtain the results that the Sivers distribution func- 
tion has the opposite signs with the factor 3 difference in magnitude for the two models, whereas 
the Boer-Mulders distribution function has the same sign and magnitude. We can understand 
these results from the properties of the light-cone wavefunction representations of the Sivers and 
Boer-Mulders distribution functions. 
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1 Introduction 



It was found that the final-state interaction of quark and gluon induces the single-spin asymmetry 
in the semi-inclusive deep inelastic scattering at the twist-two level [1]. Then, this time-odd twist- 
two effect was interpreted as the Sivers effect [2] by finding that the final-state interaction can 
be treated as the source of the time-odd Sivers distribution function [3, 4, 5, 6, 7]. It is also 
often referred to as "naively T-odd", because the appearance of this function does not imply a 
violation of time-reversal invariance, since they can arise through the final-state interactions. With 
these developments, the existence of the Sivers distribution function has gained a firm theoretical 
support. The Sivers distribution function describes the difference between the momentum 
distributions of unpolarized quark inside the nucleons transversely polarized in opposite directions. 
There is another quark distribution function of the nucleon induced by the final-state interaction 
of quark and gluon, which is called the Boer- Mulders distribution function hi [8]. h± describes the 
difference between the momentum distributions of the quarks transversely polarized in opposite 
directions inside unpolarized nucleon. The distribution functions and are depicted in Figs. 
1 and 2. 




Figure 1: Schematic depiction of the Sivers distribution function f^. The spin vector St of the 
nucleon points out of and into the page, respectively, and fcy is the transverse momentum of the 
extracted quark. 




Figure 2: Schematic depiction of the Boer-Mulders distribution function hj . The spin vector St 
of the quark points out of and into the page, respectively, and Ut is the transverse momentum of 
the extracted quark. 

The light-cone wavefunctions are valuable for studying the hadronic processes by treating the 
non-perturbative effects in a relativistically covariant way [9, 10]. The formulas which express 
the electromagnetic form factors and the generalized parton distribution functions in terms of the 
light-cone wavefunctions were found in Refs. [11, 12, 13] and Refs. [14, 15], respectively. In Ref. 
[16] the light-cone wavefunction representation of the nucleon electric dipole moment was found 
by introducing the complex phases of the light-cone wavefunctions, and studied a general relation 
connecting nucleon electric dipole and anomalous magnetic moments. 

In this paper we derive the formulas which express the Sivers and Boer-Mulders distribution 
functions in terms of the matrix elements of the nucleon spin states. Then we derive the light-cone 
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wavefunction representations of the Sivers and Boer-Mulders distribution functions, and we calcu- 
late these functions for the scalar and axial-vector diquark models by using these representations. 

The light-cone wavefunction representations show that the Sivers distribution function is given 
by the overlap of the wavefunctions of the same quark spin states of the opposite nucleon spin 
states, whereas the Boer-Mulders distribution function is given by the overlap of the wavefunctions 
of the opposite quark spin states within a given nucleon spin state. From these properties of their 
light-cone wavefunction representations, we can understand why the Sivers distribution function 
has the opposite signs with the factor 3 difference in magnitude for the scalar and axial-vector 
diquark models, whereas the Boer-Mulders distribution function has the same sign and magnitude 
for these diquark models. 



2 Sivers and Boer-Mulders Distribution Functions 

The /^-dependent unpolarized quark distribution function fi(x, k±), the Sivers distribution func- 
tion fjj,(x,k±) and the Boer-Mulders distribution function h±(x,k±) are parts of the proton 
correlation function $(rr, k± : P, S) [8]: 



$0, k ± : P, S) 



M 
2P+ 



fi(x,k A 



7 ■ P 



/11/pcr 



~W 2 



+ hi(x, k±) 



g auk i P u 



M 2 



(1) 

from which we find that fi(x,k±) and f^(x,k±) can be defined through matrix elements of the 
bilinear vector current: 



d y d 2 V± e ixP+ 

167T 3 



y--^ (p, sj^(o) 7 + Hv) \p, s±) 



y+=0 



(2) 



2P~ 
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(3) 



For an explicit calculation, let us consider the case of S± = (S]_, Sj_) = (0, 1) for the transverse 
spin in (2). Then, the proton state is given by (|P,t) +i\P, l))/y/2 and Eq. (2) becomes 



B 

where 



(P,t|-i<*UI T 



V2 



m i + ^(y) 



\P,t) + i\P,±) 



V2 



y+=0 



k 1 

fi(x,k±) - S 2 ± -j f^(x,k 



B = 



d y d V± e ixP+y--ik ± -y ± 
167T 3 



(4) 



(5) 



From (4) we have 



k]_ L 



B 



2 L 



(P,t \J + (y)\PA) + (P,i\J + (y)\P,i) 
(P,t \J + (y)\P,i) - (P,i\J + (y)\PA) 



y+=0 



y+=0 



(6) 
(7) 



where J + (y) = ip(0) ^ + ip{y). 
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On the other hand, from (1) the Boer- Mulders distribution function hi(x, k±) can be defined 
through matrix elements of the bilinear tensor current: 



dy d 2 V± e ixP+ 
16vr 3 



v~-*^ (P,S x \i,(0)a i+ iP{y) \P,S ± ) 



y+=0 



2P+ 



hi(xM)j±U(P,S ± ) 7 + U{Pj s 



which gives 
k 



M 



K (x, k±) 



B 



y+=o 



• (9) 



3 Light-Cone Wavefunction Representations of Sivers and 
Boer-Mulders Functions 

The expansion of the proton eigensolution | xjj p ) on the eigenstates { | n) } of the free Hamiltonian 
Hlc gives the light-cone Fock expansion: 



En 

n i=l 



dxi d 2 k±i 
'xll6n 3 ' 



(10) 



i=l 



i=l 



x ij n (xi, k ±i , Xi) 



n\ XiP~^~ , X{P± + k±i, A, 



The plus component momentum fractions Xj = kf / P + and the transverse momenta kj_i of par- 
tons represent the relative momentum coordinates of the light-cone wavefunctions. The physical 
transverse momenta of partons are p±i = XiP± + k±i. The Xi label the light-cone spin projections 
of the partons along the quantization direction z. The n-particle states are normalized as 



P^,p'u, K h vt.V^ \) = II 1Q ^Ptm + - P i + ) S {2) (p'± t - Pu) ■ ( n ) 



i=l 



From (6) and (7) we get 

fi{x,k L )=C~ i/jJ n *(xi,k ±i ,Xi) ipl n) {xi,k ±i ,\i) + ^{xi,k ±i ,\i) rf n) {xi,k ±i ,\i) , (12) 

k 1 - i r - -> - - 1 

-jjfi T (x,k ± ) = C- 4>l n *(xi,k_Li,Xi) rf n )(xi,k_Li, Xi) - ^(xi,k ±i ,Xi) i;j n) (xi,k±i,Xi) , (13) 

where 

c - e /n^r^ i***s( * (2) (x>^ s(x- xi )s^(k ± -k ±1 ). (M) 



n,\i i=l 



As we see in (13), the Sivers distribution function is given by the product of the light-cone 
wavefunctions which have opposite proton spin states and same quark spin states. 
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From (9) we have 



M 



(x, kj 



C 



, (-«) ( 1>l n )(?i,kj_i, Ai =|, A^i) ^f n) (:ri,A;jLi,Ai =t, A^i) 

- ^J n J(^i, fc±i, Ai =|, A^i) ^f n) (xi, fc^, Ai =|, A^i) 

+ i>(n)( X i> kj -^ ^1 ^0 ^(n)( X ii fc ±i, = t' Aj^l) 

- ^( n ) ( z i> Ai =|, A^i) ^ n) (x 4 , fc±i, Ai =|, A^i) 



(15) 



As we see in (15), the Boer-Mulders distribution function is given by the product of the light- 
cone wavefunctions which have same proton spin states and opposite quark spin states, whereas 
we found in (13) that the Sivers distribution function is given by the product of the light-cone 
wavefunctions which have opposite proton spin states and same quark spin states. 



4 Explicit Calculations in Diquark Models 
4.1 Scalar Diquark Model 

In this subsection we calculate the Sivers and Boer-Mulders distribution functions of the scalar 
diquark model by using the light-cone wavefunction representations derived in section 3. In the 
scalar diquark model, the J z = +| two particle Fock state is given by [1, 13] 



two particle 



(P+=1,Pl = 0_ 



x(l — x) I67H 



+ x, k ± ) +i[f\(x,jcj 



x , kj 



(16) 



where 



^ +l (x,k ± ) = ^l ¥ > [ 



(17) 



The scalar part of the wavefunction tp is given by 

9 1 



<p(x, k A 



x 



l-x 



k'] +B 



where 



B = x(l — x] 



■M 2 



m 



X 2 



x 1 — x 

Similarly, the J z = — \ two particle Fock state is given by 



(19) 



two particle 



(P+ = 1,P 

dxd 2 k 1 







a/x(1 — X) 167T 3 



Z 2 



2 i x 1 



(20) 
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where 

^l(xj ± ) = &2&<p. (21) 

~ 2 

The coefficients of <£> in Eqs. (17) and (21) are the matrix elements of u( ~ k+ J^ k± - U ^ P+ J^ ^- which 
are the numerators of the wavefunctions corresponding to each constituent spin s z configuration. 




Figure 3: (a) Tree level diagram and (b) diagram with final-state interaction. 



In Ref. [1] it was found that the contributing amplitudes for j*P ~> q(qq)o are given by the 
following formulas through one loop order which is depicted in Fig. 3: 

A(f^) = ^Ml C (h + i e £- gi ) (22) 
Ato-H) = C(h + &hl, (23) 

X 57T 

A(^t) = ~ t^+m C {h + i e -^g 2 ) (24) 

x Sir 

j / 1 , ,v (m + xM) .eie 2 , 

-4(^-4) = " C + ^i) , 25 



where 



and 



C = -gei P + Vx2x(l-x) (26) 

1 

fci+x(l-x)(^M 2 + ^ + £;) 



h = : ~ — . ^ 2 . A2 ; , (27) 



1 -i 

9i = I da - , (28) 

a(l-a)ki + a\ 2 g + (l-a)B 

" x -a 

g 2 = I da - . (29) 

o a(l - a)k\ + a\ 2 + (1 - a)B 
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In the above, e\ and e 2 are the quark and diquark charge, and M, m, X and X g are the nucleon, 
quark, diquark and gluon mass, respectively. Ref. [1] fixed = — Cpas, where Cp = | in order 
to relate the above results to QCD. We take X g = at the end of the calculation. We note that 
the results (28) and (29) are for the semi-inclusive deep inelastic scattering, and the results for 
the Drell-Yan process have opposite signs compared to (28) and (29) [3, 17]. 

The final-state interactions in semi-inclusive deep inelastic scattering are commonly treated as 
a part of the proton distribution function [3, 6]. If we adopt the same treatment for the wave- 
functions, we can consider that the final-state interactions for the scalar diquark model depicted 
in Fig. 3 induce the spin-dependent complex phases to the wavefunctions [18]: 



^J(x,^) = ^^(l + m 1 )^, 



(30) 



(31) 



where a± and a 2 are given by 



a h2 = ^(k 2 ± + B) g 1>2 (32) 



with <7i ;2 given in (28) and (29). 

Using the wavefunctions (30) and (31) in the formulas (12), (13) and (15), we obtain 



/ l( .,S0 = J^[(M+^ + §],*, (33) 

The results in (34) and (35) agree with those in Refs. [1, 6] with an additional overall minus sign 
which should be corrected [19]. 



4.2 Axial- Vector Diquark Model 

Jakob et al. [20] studied the scalar (s) and axial- vector (a) diquark models using the following 
nucleon-quark-diquark vertices: 

T s = g s (k 2 ), T ^ = ^S 7 M 75, (36) 

where g s {k 2 ) and g a {k 2 ) are form factors which we take as 1 in this paper for simplicity. We can 
then obtain the light-cone wavefunctions of scalar and axial-vector diquark models from Fig. 4. 
In Ref. [20], 75 (7^ + ^) appears at the vertex of the axial- vector diquark model instead of just 
7^75 appearing in (36). However, (36) is equivalent to the vertex of Ref. [20] since 4f vanishes 
when the polarization sum 

Ea4 (A) ^ A) = — 9nu + P^Pu/M 2 , which Ref. [20] used, is multiplied. 
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Figure 4: Diagram giving the light-cone wave] 'unctions of scalar and axial-vector diquark models. 



In order to obtain the light-cone wavefunctions of axial-vector diquark model, we decompose 
the polarization sum ej* = —g^ + P^Pv/M 2 of Ref. [20] to the following three polarization 
vectors e^ x ': 



( 



r 0(+l) ,1(+1) -2(+l) ,3(+l^ 



v/2 



(0,-l,-z,0), 



-^(0,+l,-z,0), 



(37) 



p3 pO 

^ (0) = fe.M,^). 



M 



M- 



In order to calculate the effects of the final-state interactions, we use for the gauge-field coupling to 
the axial- vector diquark in Fig. 3 the simple form ie2 g a/3 ((P — I) + (P — k))^ , which is equivalent, 
for each polarization state, to the gauge-field coupling to a scalar diquark [21]. We motivate this 
simple coupling by assuming that the QCD coupling to the diquark is independent of the spin 
state of the diquark. 

The two-particle Fock state for proton with J z = +1 (positive helicity) has six possible spin 
combinations for the quark and axial- vector diquark: 



two particle 



iclc (P+,P ± = 

^tl +1 (z,£±) 



dxd 2 k\ 



2 



+ ^_i +1 (x,k ± ) 



J a/x(1 — x) 167T 3 

+i +1; xP + , k ± ^ +il>li_ 1 (x,k ± 
+ \ 0; xP + , k L ^ +ip\_ l {x,k±.) 
+- -1; xP + , k L \ +^_i_ 1 (x,k 1 _ 



(38) 



-- +1; xP + , k A 



\ 0; xP~ 



1 ; xP + , k 

2 



where the two-particle states |sf , sg; x,k±) are normalized as in (11). sf and denote the z- 
component of the spins of the constituent fermion and boson, respectively, and the variables x 
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and k± refer to the momentum of the fermion. The wavefunctions are given by [21] 



Ax,k±) 

2 ~ r± 

4>ii_ n {x,k±) - 



ip\_Ax,k±) 



+ 



V'L_ 1 (a;,fc±) 








2 (-fci+ifc 2 ) 

3 a; 



za 2 



1 (m+iM) 
3 x 



1 (+fc 1 +ifc 2 ) 



1 + z'ai ) ip , 
1 + m 2 ) y? 



(39) 



where the scalar part of the wavefunction ip(x, k±) 
Similarly, the wavefunctions for proton with J z 



is given in (18). 

= — | (negative helicity) are given by [21] 



iptx +1 (x,kj 
x, k i 








Via 



(V 



2 u 

^iiAx,k ± ) 



1 (-fci+rfc 2 ) 
3 x 

1 (m+xM) 
3 a; 

2 (m+xM) 

3 a; 



1 + za 2 J , 
1 + iaA if> , 
1 + mi J if , 

+V / i™ , ( 1 + ! °^- 



(40) 



In Eqs. (39) and (40) a\ and a 2 are given by (32), and the coefficients of if are the matrix elements 
7 • e* U ^ P )p,' P± ' > which are the numerators of the wavefunctions corresponding to each 



of 



u(k + ,k ,k±) 



Vk+ 1 VP+ 

constituent spin s z configuration. 

Using the wavefunctions given in (39) and (40) in the formulas (12) and (13), we obtain 



fi(x, k A 



16tt 3 



(M + 



x 



k\ 

X 2 J 



fi T (x, k±) 



1 1 



3 16tt 3 



M 

2 _ (M+-) if 



x 



x- 8vr 1 1 ' k 2 



(fcj + B) 
B 



(41) 



(42) 



In the same way, by using the wavefunctions given in (39) and (40) in the formula (15), we obtain 
the Boer-Mulders distribution function as 



167T 3 



2 — (M+-)<f 



e±e 2 



x 



X 



(kl + B)-\n 

n k 2 



1 ,(*i + 5) 
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(43) 



The Sivers distribution function f^ T is given by the overlaps of the proton wavefunctions of 
positive and negative helicities [1, 6]. As we can see in (39) and (40), only the wavefunctions with 

= contribute to the overlaps for f^. We find that the wavefunctions with s£ = in (39) 
and (40) have the exactly same structures as the wavefunctions for the quark and scalar diquark 
system given in (30) and (31), except for the overall constant factor — ^= for the positive helicity 

wavefunctions and for the negative helicity ones. These constant factors can be understood 





Sivers Function 


Boer-Mulders Function 


Scalar Diqurrk Model 


-1 


-1 


Axial- Vector Diquark Model 


+ 1 


-1 



Table 1: Relative signs and magnitudes of Sivers and Boer-Mulders distribution functions in scalar 
and axial-vector diquark models. 

by the Clebsch-Gordan coefficients for the combination of spin | and spin 1 states [21]. Therefore, 
the overlaps of positive and negative helicities for the quark and axial- vector diquark system have 
an overall factor of — | compared to the overlaps of (30) and (31) for the quark and scalar diquark 
system. This is the reason why, in the diquark models with which we work in this paper, we have 
the same Sivers distribution functions for the cases of the scalar diquark and axial- vector diquark, 
except for the difference of the additional overall constant factor — ~ of the axial- vector diquark 
model [21], as we can see in (34) and (42). 

On the other hand, the Boer-Mulders distribution function hj^ is given by the overlaps of the 
wavefunctions of opposite quark spin states within a given proton helicity state. Concerning the 
two quark spin states for a fixed value within a given proton helicity state, the wavefunctions 
in (39) and (40) have the same structures as the wavefunctions of the scalar diquark model given 
in (30) and (31). This structure is given by the relation between the light-cone and Bjorken-Drell 
spinors [21]. Therefore, the Boer-Mulders distribution of the axial- vector diquark model is the 
same as that of the scalar diquark model as we can see in (35) and (43). We summarize the 
results in Table 1. We note that Ref. [22] studied the Boer-Mulders distribution functions 
for a variety of phenomenological models and found that the signs of are all negative for the 
models which they studied. 



5 Conclusion 

In this paper we find the light-cone wavefunction representations of the Sivers and Boer-Mulders 
distribution functions. A necessary condition for the existence of these representations is that the 
light-cone wavefunctions have complex phases. We induce the complex phases by incorporating 
the final-state interactions into the light-cone wavefunctions in the scalar and axial-vector diquark 
models, and then we calculate explicitly the Sivers and Boer-Mulders distribution functions from 
the light-cone wavefunction representations. The results are the same as those obtained from the 
direct calculation of the hadronic tensor without employing the concept of the light-cone wavefunc- 
tion, since the essential interpretation of the final-state interaction is identical in both calculations. 
However, the analysis in this paper by using the light-cone wavefunction representations is useful 
for understanding the natures of the Sivers and Boer-Mulders distribution functions in a system- 
atic way. In the light-cone wavefunction representations the Sivers distribution function is given 
by the overlap of the wavefunctions of the same quark spin states of the opposite nucleon spin 
states, whereas the Boer- Mulders distribution function is given by the overlap of the wavefunctions 
of the opposite quark spin states within a given nucleon spin state. From these properties of the 
light-cone wavefunction representations, we can understand why the Sivers distribution function 
has the opposite signs with the factor 3 difference in magnitude for the scalar and axial-vector 
diquark models, whereas the Boer-Mulders distribution function has the same sign and magnitude 
for these diquark models. 
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